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We review the theory of loss networks, including recent results on their dynamical 
behaviour. We give also some new results. 
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1 Introduction 

In a loss network calls, or customers, of various types are accepted for service provided that 
this can commence immediately; otherwise they are rejected. An accepted call remains 
in the network for some holding time, which is generally independent of the state of the 
network, and throughout this time requires capacity simultaneously from various network 
resources. 

The loss model was first introduced by Erlang as a model for the behaviour of just a 
single telephone link (see Brockmeyer et al. (1948). The typical example remains that of 
a communications network, in which the resources correspond to the links in the network, 
and a call of any type requires, for the duration of its holding time, a fixed allocation 
of capacity from each link over which it is routed (Kelly, 1986). This is the case for 
a traditional circuit-switched telephone network, but the model is also appropriate to 
modern computer communications networks which support streaming applications with 
minimum bandwidth requirements (Kelly et al., 2000). There are also other examples: for 
instance, in a cellular mobile network similar capacity constraints arise from the need to 
avoid interference (Abdalla and Boucherie, 2002). 

The mathematics of such networks has been widely studied, with interest in both equi- 
librium and, more recently, dynamical behaviour. Of particular importance are questions 
of call acceptance and capacity allocation (for example, routing), with the aim of ensur- 
ing good network performance which is additionally robust with respect to variations in 
network parameters. Call arrival rates, in particular, may fiuctuate greatly. An excellent 
review of the state-of-the-art at the time of its publication is given by Kelly (1991) — see 
also the many papers cited therein, and the later survey by Ross (1995). 

We take as our model the following. Let TZ denote the finite set of possible call, or 
customer, types. Calls of each type r ^ TZ arrive at the network as a Poisson process 
with rate v^, and each such call, if accepted by the network (see below), remains in it for 
a holding time which is exponentially distributed with mean /i^^. We shall discuss later 
the extent to which these assumptions, in particular the latter, are necessary. Calls which 
are rejected do not retry and are simply considered lost. All arrival processes and holding 
times are independent of one another. We denote the state of the network at time t by 
n{t) = {nr{t), r ^ TZ), where nr{t) is the number of calls of each type r in progress at 
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that time. The process n(-) is thus Markov. It takes values in some state space TV C Z;^, 
where R = \TZ\. We assume J\f to be defined by a number of resource constraints 

"^Ajrnr<Cj, (1) 

indexed in a finite set JT, where the Ajr and the Cj are nonnegative integers. Typically 
we think of a call of each type r as having a simultaneous requirement, for the duration 
of its holding time, for Ajr units of the capacity Cj of each resource j; however, we show 
below that the resource constraints ^ can also arise in other ways. As noted above, 
in applications of this model to communications networks, the network resources usually 
correspond to the links in the network, and when discussing the model in that context we 
shall generally find it convenient to use this terminology. We shall also find it helpful to 
define, for each r ^ TZ, the parameter ; many quantities of interest depend on 

and fir only through their ratio Kr- 

We shall say that a network is uncontrolled whenever calls are accepted subject only to 
the condition that the resulting state of the network belongs to the set J\f. Uncontrolled 
networks are particularly amenable to mathematical analysis and are in certain senses very 
well-behaved. In addition, such a network has the important insensitivity property: the 
stationary distribution of the process n(-) is unaffected by the relaxation of the assumption 
that the call holding time distributions are exponential, and depends on these holding time 
distributions only through their means. This is essentially a consequence of the detailed 
balance property considered in Section 12.11 

However, as we shall also see, the performance of uncontrolled networks may be far from 
optimal. A more general control strategy is given by requiring that a call of type r, which 
arrives when the state of the network (immediately prior to its arrival) is n, is accepted 
if and only if n G Ar for some acceptance set Ar- The sets Ar may be chosen so as to 
optimise, in some appropriate sense, the network's performance. Such networks do not in 
general possess the insensitivity property described above. 

Of interest in a loss network are both the stationary distribution it and the dynamics of 
the process n(-). For the former it is usual to compute, for each r, the stationary blocking 
probability Br, that a call of type r is rejected; here we shall find it slightly easier to work 
with the stationary acceptance (or passing) probability Pr = 1 — Br- We note immediately 
that, by Little's Theorem, the stationary expected number of calls of each type r in the 
network is given by 

where defined above. Thus acceptance probabilities may be regarded as one of 

the key performance measures in the stationary regime. 

It will be convenient to refer to the above model of a loss network — in which arriving calls 
have fixed resource requirements and in which the only control in the network is the ability 
to reject calls — as the canonical model. When considering communications networks, it 
is natural to extend this model by allowing also the possibility of alternative routing, in 
which calls choose their route according to the current state of the network. Here the 
state space should properly be expanded to record the number of calls of each type on 
each route (but see below). We consider such models in Section [3T2l 

In the case where we allow not only alternative routing, but also repacking of calls already 
in the network, the model simplifies again, and it is once more only necessary for the state 



2 



space to record the number of calls of each type in progress. Consider the simple example 
of a communications network consisting of three links with capacities C(, C2, C3, and 
three call types, in which calls of each type r = 1,2,3 require either one unit of capacity 
from the corresponding link r or one unit of capacity from each of the other two links (in 
each case the distribution of the call holding time is assumed to be the same) . If repacking 
is allowed, the state of the system may be given by n = (ni, 712,713) as usual, and it is 
easy to check that a call of any type may be accepted if and only if the resulting state of 
the network satisfies the constraints 

rir + rir' < C'j. + C'^i, r ^ r', r, r' G TZ. 

This is therefore an instance of the uncontrolled network discussed above, in which the 
coefficients Ajr and Cj must be appropriately defined. 

Exact calculations for large loss networks typically exceed the capabilities of even large 
computers, and we are thus led to consider approximations. Mathematical justification 
for these approximations is usually based on asymptotic results for one of two limiting 
schemes. In the first, which we shall refer to as the Kelly limiting scheme (see Kelly, 
1986), the sets TZ, J ^ the matrix A = (Ajv)) and the parameters /i^ held fixed, while 
the arrival rates Vr and the capacities Cj are all allowed to increase in proportion to a 
scale "parameter N which tends to infinity. In the second, which is known as the diverse 
routing limit (see Whitt, 1985, and Ziedins and Kelly, 1989), the capacity of each resource 
is held constant, while the sets TZ and J" (and correspondingly the size of the matrix A) 
are allowed to increase, and the arrival rates for call types requiring capacity at more than 
one resource to decrease, in such a way that the total traffic offered to each resource is 
also held constant (in particular, this requires that the arrival rate for any call type that 
requires capacity at more than one resource becomes negligible in the limit). Results for 
the latter scheme in particular are used to justify assumptions of independence in many 
approximations. 

For each time t, define m{t) = {mj{t), j G J), where mj{t) denotes the current occu- 
pancy, or usage, of each resource j in a loss network. Define also n' to be the stationary 
distribution of the process m(-). In particular, for the canonical model defined above, for 
each t, 

= XI ^jrnr{t); (3) 

ren 

here the process m(-) takes values in the set 

M = {m€Z-l:0<mj< Cj, j G J}, (4) 
where we write J = \J\, and the distribution tt' is given by 

7r'(m) = ^2 ■""C*^)) ^ (5) 

n: An=m 

In general the process m(-) takes values in a space of significantly lower dimension than 
that of the process n{-). This is especially so in models of communications networks which 
incorporate alternative routing. It is a recurrent theme in the study of loss networks 
that, in general, at least approximately optimal control of a network is obtained by basing 
admission decisions and, in communications networks, routing decisions, solely on the state 
of the process m(-) at the arrival time of each call. Further, in this case, a knowledge of 
the distribution tt' is sufficient to determine call acceptance probabilities. We shall also see 
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that good estimates of tt' are generally given by assuming its (approximate) factorisation 
as 

7r'(m) = n^^K)' (6) 

where each vr^- is normalised to be a probability distribution. This is a further recurrent 
theme in the study of loss networks. 

In Section [2] we consider the stationary behaviour of uncontrolled networks, reviewing 
both exact results and approximations for large networks. Our approach is based on the 
use of an elegant recursion due to Kaufman (1981) and to Dziong and Roberts (1987) 
which delivers all the classical results in regard to, for example, stationary acceptance 
probabilities, with a certain simplicity. 

More general networks are studied in Section [3l In Section 13.11 we study the problem 
of optimal control in a single-resource network, where a reasonably tractable analysis of 
stationary behaviour is again possible, and where we show that either exactly or approx- 
imately optimal control may be obtained with the use of strategies based on reservation 
parameters. In Section [3.21 we consider multiple-resource networks, allowing in particular 
the possibility of alternative routing. We again derive approximations which are known 
to work extremely well in practice. In Section S] we consider the dynamical behaviour of 
large loss networks. This is important for the study of the long-run, and hence also the 
equilibrium, behaviour of networks in the case where a direct equilibrium analysis is im- 
possible. The study of network dynamics is also the key to understanding their stability. 
Finally, in Section [5] we mention some wider models and discuss some open problems. 

2 Uncontrolled loss networks: stationary behaviour 

We study here the stationary behaviour of the uncontrolled network introduced above, in 
which calls of any type are accepted subject only to the condition that the resulting state n 
of the network belongs to the state space M defined by the capacity constraints ([1]). In 
particular we shall see, in Section [2.31 and subsequently, that most quantities of interest, in 
particular acceptance probabilities, may be calculated, exactly or approximately, without 
the need to calculate the full stationary distribution vr of the process n(-). 

2.1 The stationary distribution 

For each r & TZ, let Sr be the vector whose rth component is 1 and whose other components 
are 0. Recall that, under the assumptions introduced above, n(-) is a Markov process. For 
n,n — 5r € and r & TZ, its transition rates between n and n — dr are rirfJ-r and f^. It 
thus follows that the stationary distribution tt of the process n(-) is given by the solution 
of the detailed balance equations 

TT{n)nrfJ.r = TT{n — Sr)i^r, r G TZ, n^Af, (7) 

where, here and elsewhere, we make the obvious convention that 7r(n — Sr) = whenever 
Ur = 0. That is, 

vr(n) = G-in^' ^eAA, (8) 

r€TZ 

where the normalising constant is determined by the requirement that X^tigA^ '^('^) ~ 
1. The simple product form of the stationary distribution dH]) is a consequence of the fact 
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that the equations ^ do have a solution, that is, it is a consequence of the reversibihty of 
the stationary version of the process n(-). Note also that here the stationary distribution vr 
depends on the parameters and fir only through their ratios Kr = lyr/fJ-r, r ^ TZ. This 
result is not in general true for networks with controls. 

In the variation of our model in which calls of each type r have holding times which are 
no longer necessarily exponential (but with unchanged mean n^^), it is well-known that 
the stationary distribution vr of the process n(-) continues to satisfy the detailed balance 
equations ^ and hence also ([8]) . For a proof of this insensitivity property, see Burman et 
al. (1984). 

The stationary probability that a call of type r is accepted, is given by 

Pr= Yl ^(^)' (9) 

where Mr = {n ^ J\f : n -\- 6r E N}. In Section 12.31 we give a recursion which permits a 
reasonably efficient calculation of the probabilities Pr in networks of small to moderate 
size. However, the exact calculation of acceptance probabilities is usually difficult or 
impossible in large networks. We shall therefore also discuss various approximations. 



2.2 The single resource case 

Consider first the case TZ = {1} of a single call type. For convenience we drop unnecessary 
subscripts denoting dependence on r G 7^; in particular we write k = vj fi. We then have 
M = \n: n < C\ for some positive integer C The stationary distribution vr is a truncated 
Poisson distribution, and the stationary acceptance probability P is given by Erlang's 
well-known formula, that is, by P = 1 — 7r(C) = 1 — £"(«, C), where 

c) = ' , - (10) 

Note also that, from ([2]), the expected number of calls in progress under the stationary 
distribution vr is given by kP . 

While exact calculation of blocking probabilities via Erlang's formula (jlOp is straightfor- 
ward, it nevertheless provides insight to give approximations for networks in which C and 
K are both large. Formally, we consider the Kelly limiting scheme in which C and k are 
allowed to tend to infinity in proportion to a scale parameter A'^ with p = C / n held fixed. 
The cases p > 1, p = \ and p < 1 correspond to the network being, in an obvious sense, 
underloaded, critically loaded, and overloaded respectively. A relatively straightforward 
analysis of ([TO]) shows that, 

P min(l,p) as ^ oo. (11) 

For p > 1 the error in the approximation P ~ 1 may be estimated by replacing the 
truncated Poisson distribution of n by a truncated normal distribution: for p > 1 it may 
be shown to decay at least exponentially fast in A^, while for the critically loaded case 
p = 1 it may be shown to be 0{N~'^/'^) as N ^ oo. For the overloaded case p < 1 the 
approximation P ^ p may be refined as follows. Observe that in this case, and since k 
and C are large, it follows from either ([7]) or ([8]) that the stationary distribution of free 
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capacity in the network is approximately geometric and so the stationary expected free 
capacity is given by the approximation 

C-B^in)^-^. (12) 
1-p 

Combining this with ^ leads to the very much more refined approximation for the sta- 
tionary acceptance probability given by 

P^P-—^. (13) 

The error in this approximation may be shown to be o{N~^) as — > oo, so also that in 
the original approximation P ~ p is 0{N~^). 



2.3 The Kaufman-Dziong-Roberts (KDR) recursion 

For the general model of an uncontrolled network, we now take the set to be given by 
a set of capacity constraints of the form ([1]). We give here an efficient recursion for the 
determination of stationary acceptance probabilities, due in the case ^ = {1} to Kaufman 
(1981) and in the general case to Dziong and Roberts (1987). 

Recall that vr' is the stationary distribution of the process m{-) defined in the Introduction. 
Since a call of type r arriving at time t is accepted if and only if mj{t—) + Ajr < Cj for 
all j such that Ajr > 1 (where m{t—) denotes the state of the process m(-) immediately 
prior to the arrival of the call) , it follows that a knowledge of vr' is sufficient to determine 
stationary acceptance probabilities. Typically the size J of the set J' is smaller than the 
size R of the set TZ, and so the dimension of the space M defined by ^ is smaller than 
that of A^. Thus a direct calculation of vr', avoiding that of vr, is usually much more efficient 
for determining acceptance probabilities. 

For each r € 7^, define the vector Ar = {Ajr, j G J). For each m £ M. and r G TZ, 
summing the detailed balance equations ([7]) over n such that An = m and using also ([5]) 
yields 

KrTr'{m - Ar) = B{nr I m)7r'(m), r e TZ, m G M, (14) 

where 

tj[nr I m) — -= -— 

is the stationary expected value of rir given An = m. Since, for each m and each j, we 
have X^rGT?. I ~ "^i' follows from (I14p that 

AjrKrir'^m — Ar) = mjvr'(m), m G M, j G J. (15) 

This is the Kaufman-Dziong-Roberts (KDR) recursion on the set A4, enabling the direct 
determination of successive values of vr'(m) as multiples of vr'(O). The entire distribution 
vr' is then determined uniquely by the requirement that ^^^^j^ vr'(m) = 1. 



2.4 Approximations for large networks 

We now suppose that k^, r G 7^, and Cj, j G are sufficiently large that the exact calcula- 
tion of the stationary distributions vr or vr' is impracticable. We seek good approximations 
for the latter and for acceptance probabilities. 
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A simple approximation We give first a simple approximation, due to Kelly (1986), 
which generalises the approximation P ~ min(l,p) of Section [2.21 for the single-resource 
case. To provide asymptotic justification we again consider the Kelly limiting scheme, in 
which the parameters Kr and Cj are allowed to increase in proportion to a scale param- 
eter N, the sets TZ, J and the matrix A being held fixed. We assume (this is largely for 
simplicity) that the matrix A is such that for each m € there is at least one n G Z 
such that An = m. (This implies in particular that the matrix A is of full rank.) We 
outline an argument based on the equations (jl4p and the KDR recursion (jlSp . 

Suppose that 7r'(m) is maximised at m* € hA. The distribution dH]) of vr is a truncation 
of a product of independent Poisson distributions each of which has a standard deviation 
which is 0(A^^/^) as the scale parameter N increases. From this and from the mapping 
of TT to vr', it follows that all but an arbitrarily small fraction of the distribution of vr' 
is concentrated within a region 7W* C such that the components of m G hA* differ 
from those of m* by an amount which is again increases. Further, it 

is not too difficult to show from the above condition on the matrix A that, for each r, 
E(n,. I m) varies smoothly with m, and that within M* we may make the approximation 
E(nr I m) ~ E(nr | m*) (the error yet again being 0(A^^/^) as increases). It now 
follows from (I14|) that within TM* we have 

v^'(m)«v^'(m*)^P?*""^ (16) 

where necessarily, since ra* maximises 7r'(m), 

< < 1, i G J, (17) 

■pj = 1, for j such that m* < Cj . (18) 

Further, from p5|) . 

Y.Ajr^^rUpt''' =m* <Cj, jej. (19) 

Thus, from (fT6]) . within A4* the stationary distribution n' of m(-) does indeed have 
the approximate factorisation ([6]), where each of the component distributions ttj is here 
geometric (and where in the case pj = 1 the geometric distribution becomes uniform). 
Further, for each r and each j, we have 

TTjiimj : rrij < Cj - Ajr}) ^ pf'". 

Thus the stationary acceptance probabilities Pr are given by the approximation 

Pr^llpf'\ r^TZ. (20) 

Kelly (1986) considered an optimisation problem from which it follows that the equa- 
tions (fT71l - (fT9]l determine the vectors m* and p = (pj, j G J) uniquely. He further 
showed, in an approach based on consideration of the stationary distribution vr, that the 
approximation ()20p becomes exact as the scale parameter N tends to infinity. 
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A refined approximation The (multiservice) reduced load or knapsack approximation 
(Dziong and Roberts, 1987, see also Ross, 1995) is a more refined approximation than that 
defined above. It is given by retaining the approximate factorisation ([6]) of the stationary 
distribution vr' of m(-). However, subject to this assumed factorisation, the estimation of 
the component distributions tTj is refined. 

For each j ^ J and r E 7^, define 

pjr= Yl (21) 

mj=0 

note that pjr = 1 if Ajr = 0. For fixed j, substitution of ([6]) into the KDR recursion (jlSp 
and summation over all for all k ^ j yields 

Y^ Ajri nrYl^PkrjT^jimj — Ajr) = mjTTj{mj), l<rnj<Cj, j&J (22) 

(where, as usual, we make the convention 7r^(mj) = for rrij < 0). This is the one- 
dimensional KDR recursion associated with a single resource constraint j, and is readily 
solved to determine vr^- and hence the probabilities pjr, r G 7^, in terms of the probabil- 
ities pfcD r ^ TZ, for all k ^ j. We are thus led to a set of fixed point equations in the 
probabilities pj>, for which the existence — but not always the uniqueness, see Chung and 
Ross (1993) — of a solution is guaranteed. From ([6]), the probability that a call of type r 
is accepted is then given by 

Pr=\{pjr- (23) 

We remark that the recursion (I22p corresponds to a modified network in which there is 
a single resource constraint j and each arrival rate Kr is reduced to HrY\k=/=j Pkr- This 
reduced load approximation is of course exact in the case of a single-resource network. 

In the case where each Ajr can only take the values or 1 we may set pj = pjr for r such 
that Ajr = 1. The fixed point equations (|21|) and (|22|) then reduce to 

Pj = ^-E(j2^rllpt''^Cj\ (24) 

XrelZ k^j ) 

where E is the Erlang function (llOp . This case is the well-known Erlang fixed point 
approximation (EFPA) and has a unique solution, see Kelly (1986), and also Ross (1995). 
It yields acceptance probabilities which are known to be asymptotically exact in the Kelly 
limiting scheme discussed above, and also, under appropriate conditions, in the diverse 
routing limit discussed in the Introduction — see Whitt (1985), and Ziedins and Kelly 
(1989). The EFPA also has an extension to the case of general Ajr, which may be regarded 
as a simplified version of the reduced load approximation. As with the latter approximation 
the EFPA may here have multiple solutions. 

3 Controlled loss networks: stationary behaviour 

We now study the more general version of a loss network, in which calls are subject to 
acceptance controls, and the issues are those of achieving optimal performance. 
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3.1 Single resource networks 



We consider a simple model which illustrates some ideas of optimal control — in particular 
those of robustness of the control strategy with respect to variations in arrival rates (which 
may in practice be unknown, or vary over time). 

Suppose that TZ = {1, 2} and that as usual calls of each type r arrive at rate I'r and have 
holding times which are exponentially distributed with mean fi^^- Suppose further that 
there is a single resource of capacity C and that a call of either type requires one unit of 
this capacity, so that the constraints ([T]) here reduce to ni + n2 < C. We assume that 
calls of type 1 have greater value per unit time than those of type 2, so that it is desirable 
to choose the acceptance regions 2, so as to maximise the linear function 



for some ai > 02 > (where, again as usual, for each call type r, = Vrl y^-r and is the 
stationary acceptance probability.) An upper bound for the expression in (|25p is given by 
the solution of the linear programming problem, in the variables Pi, 

maximise (/>(Pi, P2), subject to Pr G [0, 1] for r = 1, 2, kxP\ + ^2^2 < C (26) 

(where the latter constraint follows from ([2])). It is easy to see that the solution of this 
problem is characterised uniquely by the conditions 



It is clearly not possible to choose the acceptance regions A2 so that the corresponding 
values of Pi, P2 solve exactly the problem (j26p . However, we show below that this solution 
may be achieved asymptotically as the size of the system is allowed to increase, and further 
that there is an asymptotically optimal control that is both simple and robust with respect 
to variations in the parameters ki, k^- 

We consider first the form of the optimal control in the special case [i\ = [i2- Here since, 
at the arrival time t of any call, those calls already within the system are indistinguishable 
with respect to type, it is clear that the optimal decision on call admission is a function 
only of the arriving call type and of the total volume m{t—) = ni(i— ) + n2{t—) of calls 
already in the system. A formal proof is a straightforward exercise in Markov decision 
theory. Further, simple coupling arguments show that, for an incoming call of either type 
arriving at time t and any < m < C, if it is advantageous to accept the call when 
m{t—) = m, then it is also advantageous to accept the call when m{t—) = m — 1. It 
follows that the optimal acceptance regions are of the form 



for some reservation parameter k, whose optimal value depends on C, ki and K2- 
Consider now the general case where we do not necessarily have /ii = fi2, and suppose 
that C, I'l and 1^2 are large. More formally we again have in mind the Kelly limiting 
scheme in which these parameters are allowed to increase in proportion to some scale 
parameter which tends to infinity (while /ii, ^2 are held fixed). We further suppose 




(25) 



Pi = P2 = 1, whenever kiPi + K2P2 < C 
P2 = 0, whenever Pi < 1. 



(27) 
(28) 



^1 = {n: ni + n2 < C} 
A2 = {n: ni + n2 < C — k} 



(29) 
(30) 
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that the acceptance regions are again as given by ([^^ and ([HUjl . where the reservation 
parameter k increases slowly with N , i.e. in such a way that 

k^oo, k/C^O, asAT^oo. (31) 

It is convenient to let Pi, P2 denote the limiting acceptance probabilities. In the case 
K1 + K2 < C, it is not difficult to see that, since k/C ^ as N ^ 00, we have Pi = P2 = 1, 
so that Pi, P2 solve the optimisation problem (f^Hj) . Consider now the case ki + K2 > C. 
Here, again since k/C — >OasA^— > 00, it follows that, in the limit, the capacity of the 
network is fully utilised. Further, if ki is sufficiently large that Pi < 1 (informally, even 
for large N, calls of type 1 are being rejected in significant numbers), then the effect of the 
increasing reservation parameter k is such that, again in the limit, the network remains 
sufficiently close to capacity to ensure that no calls of type 2 are accepted, and hence 
P2 = 0. It now follows that when ni + K2 > C, the limiting acceptance probabilities Pi, 
P2 satisfy the conditions (j27p and (j28p and so again solve the optimisation problem (|26p . 

The above analysis demonstrates the asymptotic optimality of any strategy based on the 
use of a reservation parameter k, provided only that, in the limiting regime, k increases 
in accordance with (I3ip . In practice, in a large network (here for large C), only a small 
value of k is required in order to achieve optimal performance. We also observe that the 
performance of a reservation parameter strategy is indeed robust with respect to variations 
in Ki, K2. 

This analysis also extends easily to the case where there are more than two call types, 
and also, with a little more difficulty, to that where the capacity constraint is of the form 
J^ren^rnr < C for general positive integers (see Bean et a/., 1995). Here a different 
reservation parameter may be used for each call type, and, in the Kelly limiting scheme, a 
complete prioritisation and optimal control are again achieved asymptotically by allowing 
the differences between the reservation parameters to increase slowly. 

3.2 Multiple resource models 

Consider now the general case of the canonical model in which there is a set of resources 
and in which state n of the network is subject to the constraints ([T]). Suppose that 
it is again desirable to choose admission controls so as to maximise the linear function 
(j){P) := X]^=i ttrKrPr of the stationary acceptance probabilities Pr, for given constants a,., 
r £lZ. As in Section 13. H we may consider the linear programming problem 

R 

maximise (t>{P), subject to Pr G [0, 1] for r £ TZ, ""^^ AjrKrPr < Cj for j G J, (32) 

r=l 

which provides an upper bound on the achievable values of the objective function </>. It is 
easy to see that this value may be asymptotically achieved within the Kelly limiting regime 
by reserving capacity Ajr^rPr at each resource j solely for calls of each type r, where here 
P is the solution of the problem ()32p . However this strategy is neither optimal in networks 
of finite capacity, nor is it robust with respect to variations in the parameters Kr. At the 
opposite end of the spectrum from this complete partitioning policy is that of complete 
sharing. The latter can lead to unfairness if there are asymmetric traffic patterns, with the 
potential for some call types to receive better service than others. In practice it is expected 
that good strategies will be based on the sharing of resources and the use of reservation 
parameters — as was shown to be optimal for single resource networks in Section 13.11 
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In the case of communications networks it is natural to allow also alternative routing, as 
described in the Introduction. An upper bound for the achievable performance is given by 
supposing that repacking is possible, i.e. that calls in progress may be rerouted as necessary. 
In this case, our model for the network reduces to an instance of the canonical model (as 
defined in the Introduction) with appropriately redefined set J, matrix A = (Ajr) and 
capacities Cj. The upper bound on (p^P) given by the linear programming problem ()32p 
is then also an upper bound in the more usual case in which repacking is not allowed. 
In the latter case practical control strategies are again based on the use of appropriate 
reservation parameters, and there is some hope that performance close to the upper bound 
above may be achieved in networks with sufficiently large capacities or sufficient diversity 
of routing, even without repacking. In applications reservation parameters are generally 
used to prioritise different traffic streams. In networks with alternative routing they also 
prevent the occurrence of network instabilities, where, for fixed parameter values, the 
network may have two or more relatively stable operating regimes — one in which most 
calls are directly routed, and others in which many calls are alternatively routed, with a 
resulting severe degradation of performance (see Gibbens et al, 1990, Kelly, 1991). By 
giving priority to directly routed traffic, the use of reservation parameters prevents the 
network from slipping into an inefficient operating state. 

There have been numerous investigations of control strategies for communications networks 
that employ either fixed or, particularly, alternative routing. Such strategies are often 
studied in the context of fully connected networks. Two of the most commonly studied are 
least busy alternative (LBA) routing and dynamic alternative routing (DAR). LEA routing 
seeks to route calls directly if possible, and otherwise routes them via that path which 
minimises the maximum occupancy on any of its links. Directly routed calls are usually 
"protected" with some form of reservation parameter (Kelly, 1991, Marbukh, 1993). Hunt 
and Laws (1993) showed that, for fully connected networks which permit only two- link 
alternative routes, LBA routing is asymptotically optimal in the diverse routing limit (see 
Section . This policy is robust to changes in traffic patterns, but has the difficulty 
that it requires information on the current states of all possible alternative paths before 
an alternative routing decision is made. 

A much simpler routing scheme is DAR (Gibbens et al, 1989, Gibbens and Kelly, 1990). 
In this scheme, for each pair of nodes, a record is maintained of the current preferred 
alternative route, and this is the one that is used if a call cannot be routed directly. If 
neither the direct route nor the current preferred alternative route are available, then the 
call is rejected, and a new preferred alternative route is chosen at random from those 
available. Directly routed traffic is again usually protected by a reservation parameter. 
This policy is easy to implement. It does not require information about the current state 
of the system to be held at any node, just a record of the current preferred alternative 
route to other nodes. It is also robust to changes in traffic patterns — alternative routes 
on which the load increases will be discarded and replaced by routes on which the load 
is lower. Neither LBA routing nor DAR require traffic rates to be known or estimated 
(except approximately, in order to set the appropriate level of the reservation parameters). 

Acceptance probabilities for controlled loss networks are usually estimated using a gen- 
eralised version of the reduced load or knapsack approximation of Section 12.41 As there, 
we make the approximation ([6]) for the stationary distribution vr' of the resource occu- 
pancy process 'm{-). Each of the marginal distributions ttj is estimated as the stationary 
distribution of a Markov process on {0, . . . ,Cj} which approximates the behaviour of the 
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resource j considered in isolation. Let pjr be the probability under this distribution that 
a call of type r is accepted, subject to the controls of the model, with pjj. = 1 if Ajr = 0. 
In the case of the canonical model, in which no alternative resource usage is allowed, calls 
of each type r are assumed to arrive at resource j at a rate i^r Y\kjtj Pkr — this is the "re- 
duced load" for calls of type r at this resource; further, calls of this type arriving at this 
resource are subject to the acceptance controls of the model and, if accepted, depart at 
rate fir as usual. The estimated stationary distribution tt^- then determines the acceptance 
probabilities pjr at the resource j. Thus we are again led to a set of fixed point equations 
which determine — not always uniquely — the acceptance probabilities pjr for all r G 7?. and 
j G J . Finally the stationary network acceptance probability for calls of each type r 
is again given by Pr = Wj^^jPjr- 

In the case of a communications network where the canonical model is extended by allowing 
the possibility of alternative routing, it is necessary to modify the above approximation. 
Suppose, for example, that a link (resource) j forms part of the second choice route for 
calls of type r. Then, in the one-dimensional process associated with link j, the arrival 

rate for calls of type r is taken to be the product of the arrival rate at the network, 
the probability that a call of this type is rejected on its first-choice route, and (as before) 
the probabilities that the call can be accepted at each of the remaining resources on the 
alternative route (see e.g. Gibbens and Kelly, 1990). 

The basis of the reduced load approximation is the approximate factorisation of the dis- 
tribution tt' above. In the case of controlled networks, this approximation fails to become 
exact under the Kelly limiting regime in which capacities and arrival rates increase in 
proportion. It may, however, be expected to hold under sufficiently diverse routing. It is 
known to be remarkably accurate in most applications. 

4 Dynamical behaviour and stability 
4.1 Fluid limits for large capacity networks 

We now consider the dynamical behaviour of large networks. As well as such behaviour 
being of interest in its own right — for example in networks in which input rates change 
suddenly, fixed points of network dynamics correspond to equilibrium, or quasi-equilibrium, 
states of the network (see below). The identification of such points is often the key 
to understanding long-term behaviour, in particular to resolving stability questions and 
determining stationary distributions where (as is usual) the latter may not be directly 
calculated. However, we note that it is characteristic of loss networks that, from any 
initial state, equilibrium is effectively achieved within a very few call holding times, so that 
transient performance is of less significance than is the case for networks which permit 
queueing. 

We describe a theory first suggested by Kelly (1991). We yet again assume the Kelly 
limiting scheme described in the Introduction, in which the network topology is held fixed 
and arrival rates and capacities are allowed to increase in proportion. More explicitly, we 
consider a sequence of networks satisfying our usual Markov assumptions (though this is 
not strictly necessary) and indexed by a scale parameter N . All members of the sequence 
are identical in respect of the (finite) sets TZ, J, the matrix A = {Ajr,j ^ J, r eTZ), and 
the departure rates fir, r eTZ. For the Nth member of the sequence, calls of each type r 
arrive at rate Nur for some vector of parameters u, and the capacity of each resource j 
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is NCj for some vector of parameters C, where, for simplicity, we take each Cj to be 
integer- valued. As always, it is convenient to define Kr = Vrll^^r for each r eTZ. 
We now describe the rules whereby calls are accepted. For each N, let n^(t) = {n^{t), r G 
IZ), where {t) is the number of calls of type r in progress at time t. Define also the free 
capacity process fh^ {■) = (m^(-), j G J') where each rnj^ (t) = NCj — "^r^-ji^jr-n^ (i) is 
the free capacity of resource j at time t. A call of type r arriving at time t is accepted 
if and only if the free capacity m^(t— ) of the system, immediately prior to its arrival, 
belongs to some acceptance region Ar C Z'^. We take the acceptance regions Ar, r eTZ, 
to be independent of N, although, in a refinement of the theory, some dependence may be 
allowed. Note that, in a change from our earlier conventions, the acceptance regions Ar 
are defined in terms of the free capacity of each system. 

While the above description defines instances of the canonical model of the Introduction, 
more sophisticated controls, such as those involving the use of alternative routing in com- 
munications networks, may be modelled by the suitable redefinition of input streams and 
acceptance sets (see Hunt and Kurtz, 1994). 

For each N, define the normalised process £C^(-) = n^{-)/N, which takes values in the 
space 

X = {x eR^:Y^ AjrXr < Cj for all j e J}. (33) 

Assume that, as ^ co, the initial state x^{0) converges in distribution to some a;(0) € 
X, which, for simplicity, we take to be deterministic. Then we might expect that the 
process a;-^(-) should similarly converge in distribution to a fluid limit process a;(-) taking 
values in the space X, with dynamics given by 

Xr{t) = Xr{0) + / (fr-Pr(ii) — firXr{u))du, r E TZ, (34) 

Jo 

where, for each t, Pr{t) corresponds to the limiting rate at which calls of each type r are 
being accepted at time t. 

A rigorous convergence result is given by Hunt and Kurtz (1994). A somewhat techni- 
cal condition (always likely to be satisfied in applications) is required on the acceptance 
sets Ar- However, the main difficulty is that in some, usually rather pathological, cases 
the limiting acceptance rates Pr{t) may fail to be unique. 

In many cases, though, it is possible to show that, for each r, there does exist a unique 
function Pj- on X such that, for each t, we have Pr{t) = Pr{x(t)). In general, the trajec- 
tories of the limit process x{-) are then deterministic functions of their initial positions 
x{0). The fixed points x of the limit process x{-) are given by the solutions of 

UrPrix) = jlrXr, r E TZ. (35) 

In the case of a single fixed point x, to which all trajectories of a;(-) converge, it may be 
shown that the stationary distribution of the original normalised process a;-^(-) converges 
to that concentrated on the single point x. Then in particular, for each r, Pr{x) is the 
limiting stationary acceptance probability for calls of type r. In the case of multiple fixed 
points, those which are locally stable correspond to "quasi-stationary" distributions of the 
process x'^{-), i.e. regimes which are maintained over periods of time which are lengthy 
but finite. 
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4.2 Single resource networks 



As the simplest non-trivial application of the above theory, we consider the case J = 1 of 
a single resource, for which equilibrium behaviour was described in Section [3. II It is again 
convenient to write Ar for Air for each r, and similarly C for Ci. The technical condition 
referred to above on the acceptance sets Ar ^ Z+, here reduces to the requirement that, 
for each r, either m ^ Ar for all sufficiently large m € — we let IZ* denote the set of 
such r — or m ^ Ar for all sufficiently large m E Z^-. 

Here the functions Pr defined above always exist (see Hunt and Kurtz, 1994). To identify 
them, define, for each x € X, the Markov process rhx{') on ^+ with transition rates given 

by 

_^^{m-Ar at rate (35) 
+ at rate ^rXr, 

Let Tlx be the stationary distribution of this process where it exists. Define X C X by 

X = {x ^ X: ArXr = C and tTx exists}. (37) 

(The set X may be thought of as consisting of those points in X for which the limiting 
dynamics are "blocking".) Then, for a; G X, we have Pr{x) = T^x{Ar) for all r; for 
a; G X \ X, we have Pr{x) = 1 for r G 7^* and Pr{x) = for r ^ TZ* . The fixed points x 
of the limiting dynamics (in general there may be more than one such) are then given by 
the solutions of (1351). 



Consider now the case of reservation-type controls, and suppose that the call types are 
arranged in order of decreasing priority. The acceptance regions are thus given by Ar = 
{fh: fh > kr + Ar} for some = ki < k2 < ■ ■ ■ < and we have TZ* = TZ. It is easy to 
see that, in the light traffic case given by "^reiz ^r^r < C, the single fixed point x of the 
limiting dynamics is given by Xr = Kr for all r, and that all trajectories of these dynamics 
converge to x. In the heavy traffic case given by Ylreiz -^ri^r > C, define X C X by 

X = {x £ X : ArXr = C and Xr < Kr for all r G TZ}. 

Then it is straightforward to show that X C X and that all fixed points of the limiting 
dynamics lie within X (see Bean et al, 1995). In the case where Ar = 1 for all r, it is 
also straightforward to show that there is a unique fixed point. It is unclear whether it is 
possible, for more general Ar, to have more than one fixed point. 

Now define ro > to be the maximum value of r G 7^ such that Ylr<ro ^r^r < C. Suppose 
that the reservation parameters ki, . . . ,kr are allowed to increase. Further consideration 
of the processes tTx shows that, in the limit (formally as these reservation parameters 
tend to infinity), the fixed point x is necessarily unique and is such that Pr{x) = 1 for 
all r < ro, with, in the heavy traffic case, < Prg+i{x) < 1 and Pr{x) = for all 
1^ ^ ''0 + 2. Since the stationary distributions associated with our sequence of networks 
converge to that concentrated on the unique fixed point x, it follows that the reservation 
strategy does indeed approximate, and in the limit achieve, the complete prioritisation of 
call types discussed in Section 13.11 As mentioned there, and as easily verified from the 
above analysis, quite small values of the reservation parameters ki, . . . ,kr are sufficient to 
achieve a very good approximation to this prioritisation. 
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Even in the present single-resource case it is possible to achieve nonuniqueness of the 
fixed points of the hmit process x(-) by the use of more general, and sufficiently perverse, 
controls, in particular with the use of acceptance sets of the form Ar = {rh: < rh < 
kr + Ar} for some kr > (see Bean et al., 1997). Thus we may construct networks which 
have several (very different) regimes which are quasi-stationary in the sense discussed 
above. 



4.3 Multi-resource networks: the uncontrolled case 

We now consider multi-resource networks, and again study the behaviour of the fluid limit 
process x{-) associated with the Kelly limiting scheme. Here in general a rich variety of 
behaviour is possible. However, in the case of the uncontrolled networks of Section [21 in 
which calls of all types are accepted subject only to the availability of sufficient capacity, 
the process x{-) is rather well-behaved. Note that here, in terms of the available free 
capacity, the acceptance sets are given by, for each r G TZ, 

Ar = {rh : fhj > Ajr for all j}. (38) 
Recall also that X is as given by (I33p . Define the (real- valued) concave function f on X 

by 

f{x) = {Xr log Ur — X^ log flrXr + Xr) (39) 

and let x be the value of x which maximises f{x) in X. Kelly (1986) shows that, as 
N ^ oo, the stationary distribution of the process x^ {■) converges to that concentrated 
on the single point x. (Indeed this is the basis of his original derivation of the the limiting 
acceptance probabilities considered in Section f2.4[ ) 

Assume for the moment the unique existence of the functions Pr on X introduced above. 
Then, for the fluid limit process x{-), it follows from ()34p and ()39p that df{x{t))/dt = 
g[x[t)) where the function g on X is given by 

g{x) = ^ ^^^{vrPr{x) - flrX, 

= ^logZ^r — log HrXr^ (^l/rPrix) — fIrX, 

ren 

Analogously to the preceding section, for each x G X, the limiting acceptance probabilities 
Pr{x) are given by consideration of the stationary distribution of a "free capacity" Markov 
process whose transition rates depend on x. Some simple analysis of the equilibrium 
equations which define this stationary distribution (see Zachary, 2000) now shows that 
g{x) > for all a; G X with equality if and only if a; = a;. 

Thus the dynamics of the limit process x{-) are such that, away from the point x, the 
function f{x{-)) is always strictly increasing. It thus acts as a Lyapunov function, ensuring 
that all trajectories of the process x{-) converge to the single fixed point x. Indeed a 
rigorous application of the fluid limit theory of Hunt and Kurtz (1994) (again see Zachary, 
2000, for details) shows this result continues to hold even if the functions Pr on X are 
not uniquely defined (whether this can ever happen in the case of uncontrolled networks 
remains an open problem). The result therefore establishes an important stability property 
of uncontrolled networks, and guarantees that the stationary distribution describes the 
typical behaviour of the network. 
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4.4 Multi-resource networks: the general case 




For general multi-resource networks, the fluid limit process x(-) associated with the Kelly 
limiting scheme may fail to be unique, and may in particular exhibit multiple fixed points. 
We describe in some detail an elementary example, which is a simplification of one due to 
Hunt (1995b). Suppose that R = 3, J = 2, and that the matrix A is given by 

A = 

Thus in particular calls of types 1 and 2 each require capacity from a single resource, while 
calls of type 3 require capacity from both resources in the network. Suppose further that 
the (free capacity) acceptance sets are given by, for some ki,k2 > 1, 

Ai = {rh: 1 < rhi < ki}, A2 = {Trn: 1 < ?fi2 < k2}, Aj, = {m: fhi > 1, 7712 > 1}. 

(As Hunt remarks, this is not entirely unrealistic: in more complex networks, operating 
under some form of alternative routing, certain resources may have calls of certain types 
routed over them precisely when the network is in general very busy.) Finally suppose 
that Hr = 1 for all r and that the vectors v and C defined in Section [4. II (each to be scaled 
by N for the A^th member of the sequence of networks) are given by = (ui, 1^2, 1's) and 
C = {C,C). 

The process a;(-) takes values in the space X = {x G Mi]_ : a^i + X3 < C, X2+ < C}. Its 
dynamics may be determined through the fluid limit theory outlined above. For x € Xq := 
{x X : xi + < C, X2 + X3 < C} (corresponding to limit points of the dynamics well 
away from the capacity constraints) the limiting acceptance probabilities are well-defined 
and given by 

Pi{x) = P2{x) = 0, P3{x) = l. (40) 

For X G Xi := {x G X : xi + x-^ = C, X2 + Xs < C} and for x £ X2 := {x G X : xi + X3 < 
C, X2 + X3 = C} (corresponding in both cases to limit points of the dynamics such 
that only one capacity constraint is relevant) the limiting acceptance probabilities are 
again well-defined and given by consideration of a Markov process on Z+ as in the single 
resource case considered in Section 14.21 (For x € Xi, for example, it follows from the 
definition of A2 that the transition rates of this Markov process are as if 1^2 = 0.) For 
X € X12 := {x £ X : xi + X3 = C, X2 + X3 = C} it is necessary to consider also a "free 
capacity" Markov process on Z^. 

In the case 1^3 < C, these Markov processes all fail to possess stationary distributions 
and the limiting acceptance probabilities are given by ()40p for all a; G X. Thus the limit 
process a; (•) is as if 1^1 = 1^2 = and all trajectories of this process are deterministic 
functions of their initial values and tend to the single fixed point x = (0,0, u^). 

The case 1^3 > C is more interesting. Here it is readily verified that the limit process a;(-) 
possesses no fixed points in Xq. Within Xi consideration of the stationary distribution 
of the Markov process defined in Section 14.21 shows that there is a single fixed point 
x^^^ = (ai, 0, C — ai) for some oi which is independent of 1^2- Similarly within X2 there is 
a single fixed point x^"^^ = (0, 02, C — 02) for some 02 which is independent of ui. However, 
within X12 the dynamics of the limit process x{-) are not deterministic. It is further not 
difficult to show that all trajectories of a;(-) which avoid the set X12 tend deterministically 
to one of the two fixed points a;*^"'^'', x^'^'^ above (depending on whether the set Xi or the 
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set X2 is hit first). Those trajectories of x{-) which do hit X12 may, in an appropriate 
probabilistic sense, tend to either x^"^^ or x^'^\ 

The interpretation of the above behaviour is the following. Suppose that N is large and 
that, for example, resource 1 fills to capacity first. Then this resource remains full and 
blocks sufficient of the type 3 calls to ensure that resource 2 remains only partially utilised, 
with no calls of type 2 ever being accepted. This corresponds to a "quasi-stationary" 
state whose limit, as — > 00, is concentrated on the fixed point x^^\ Alternatively, if 
resource 2 fills to capacity first, the network settles, for an extended period of time, to a 
quasi-stationary state whose limit is concentrated on the fixed point a;*^^\ While, for finite 
N , transitions between these two quasi-stationary states will eventually occur, the time 
taken to do so can be shown to increase exponentially in N . 

The behaviour in the above example is typical of that which may occur in more general 
networks — in particular those using alternative routing strategies — which are poorly con- 
trolled. Fluid limits may be used to study behaviour in networks with high capacities 
and correspondingly high arrival rates, and to choose values of, for example, reservation 
parameters so as to ensure that the network does not spend extended periods of time in 
states in which it is operating inefficiently. A realistic example here is the fully-connected 
network with alternative routing considered in Section 13.21 

As noted above, fluid limits may also be used to study equilibrium behaviour, especially 
in the case where all trajectories of the limit process x{-) tend to a unique fixed point x. 
In particular we may show that, for the Kelly limiting regime considered here, the limiting 
stationary distribution of the free capacity processes {■) in general only has a product 
form in the case of uncontrolled networks. This product- form assumption is the basis of 
the commonly used approximations considered in Section f3.2l Its justification owes more 
to the results for the diverse routing limit also considered there and in Section [4.51 

4.5 The diverse routing limit 

In this section we consider the fluid limit obtained under the diverse routing regime dis- 
cussed in the Introduction. Although a high degree of symmetry is required in order to 
obtain formal limits, the results obtained lend support to the commonly made assumptions 
of independence of resource blocking which are used, for example, in the construction of 
the approximations discussed in Section 3.2. 

As outlined earlier, the diverse routing regime holds when the numbers of resources and 
possible "routes" in the network increase, while the total capacity and arrival rate at each 
resource remains constant. For this limit to exist we require a high degree of symmetry in 
the network. There are two canonical examples (with variants) that have been extensively 
studied. We describe both here using the terminology of communications networks. 

The first is the so-called star network (see, for instance, Whitt, 1985, Ziedins and Kelly, 
1989, Hunt, 1995a). Here there are K links, each with capacity C. The scale parameter of 
the regime is then taken lohe, K. Assume that calls of any size r > 1 require unit capacity 
at each of r resources and have holding times with unit mean. Then in a symmetric network 
there are (^) possible choices of the set of links for such a call. Let the arrival rate for each 
such choice be = Vr/i^Zi)^ that the total arrival rate at each resource for calls of 
size r is exactly A^. For example, we may assume that the K links are distributed around 
a central hub, through which all communications must pass. Many variants of this model 
are possible — multiple call sizes can coexist in the network, as can multiple capacities. 
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provided only that the proportion of hnks with any given capacity remains constant as K 
increases. The network is assumed to have fixed routing and the only permissible controls 
are those on admission. 

Let x^{t) = {x^{t), j € i7) where (t) is the proportion of links in which j units of 
capacity are in use at time t. For the network without admission controls, Whitt (1985) 
showed that, given the initial point x{0), the process x^{-) converges weakly to a deter- 
ministic limit process x{-), which satisfies a set of first-order differential equations with 
a unique fixed point x, such that x{t) ^ x as t ^ oo for all initial x{0). The limit x 
coincides exactly with that given by the Erlang fixed point approximation. Recall that 
the latter is obtained from the assumption that the stationary free capacity distributions 
on the various links of the network are independent of each other. For the case where 
all calls are of size two, Hunt (1995a) obtained a functional central limit theorem for the 
process x^{-), with the limit an Ornstein-Uhlenbeck diffusion process (as previously con- 
jectured by Whitt), which was then extended to more general sizes and initial conditions 
by Graham and Meleard (1995). In the case of networks with admission controls very little 
has been proved. MacPhee and Ziedins (1996) studied such networks and gave a weak 
convergence result for the process x^{-). However, there remain many open questions 
about the behaviour of this process. 

The second canonical example of the diverse routing regime is that of the fully connected 
network with alternative routing (Hunt and Laws, 1993). Here both admission and routing 
controls are possible. The network has N nodes; between each pair of these there is a link 
with capacity C, so that the total number of links is K = . Here again K is the scale 
parameter. Calls arrive at each link at rate u; each call has a unit capacity requirement 
and holding time of mean 1. There are three possible actions on the arrival of a call: 
(i) accept the call at that link, (ii) select a pair of links that form an alternative route 
between that pair of nodes and route the call along this, or (iii) reject the call. Hunt and 
Laws showed that an asymptotically optimal policy, in the sense of minimising the average 
number of lost calls in equilibrium, is to route a call directly if possible and otherwise to 
route it via an alternative route, provided that the remaining free capacity on each link of 
the alternative route is at least some reservation parameter k, where the optimal choice 
of k is determined by the parameters K and v. The optimal choice of alternative route 
is given by choosing that which is least busy, i.e. which maximises the minimum of the 
free capacities on the two links. The analysis of Hunt and Laws largely dispenses with 
the graph structure inherent in the choice of alternative routes, an assumption justified by 
analogy with earlier results of Crametz and Hunt (1991) in relation to the simpler model 
without reservation. 

As in the example of the star network, of interest here is the process x^{-), defined as 
earlier. Hunt and Laws showed weak convergence of this process to a deterministic limit 
process. They showed that this limit process satisfies differential equations which yield the 
constraints for a linear programming problem, the solution to which gives an upper bound 
on the acceptance probabilities. (These constraints correspond to the detailed balance 
equations that in equilibrium govern the changes in occupancy of a single link.) They 
further showed that their policy achieves this upper bound. 
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5 Further developments and open questions 



Our discussion has of necessity omitted many topics of interest, some of which we mention 
briefly here, as well as discussing some remaining open questions. 

One such topic is the application of large deviations techniques to loss networks in order 
to estimate, for example, blocking probabilities in cases where it is important to keep 
these very small. For an excellent introduction to this see, for instance, Shwartz and 
Weiss (1994); later papers include those by Simonian et al. (1997) and by Graham and 
O'Connell (2000). 

In some models of communications networks, particularly those whose graph structure is 
tree-like, the network topology may be such as to lend itself to more accurate calculations 
of acceptance probabilities, involving recursions that do not make the link independence 
assumption ([6]) that is such an essential feature of the approximations presented above 
(see Zachary and Ziedins, 1999). 

Extensions of loss network models include recent work by Antunes et al. (2005) which 
studies a variant of the model where customers may obtain service sequentially at a number 
of resources, each of which is a loss system. The aim here is to model a cellular wireless 
system where a call in progress may move from base station to base station. Several 
authors have also considered explicitly systems with time- varying arrival rates and/or 
retries (see, for example, Jennings and Massey, 1997, and Abdalla and Boucherie, 2002). 

A large number of interesting and important open problems remain. The approach to 
most of these seems to lie in a better understanding of network dynamics. There has 
been no systematic investigation of how to achieve asymptotically optimal control in a 
general network (for example in the sense of Section [3. 2p . using controls which are simple, 
decentralised, and robust with respect to variations in network parameters, although, for 
communications networks, there is a belief that this will usually combine some form of 
alternative routing with the use of reservation parameters to guarantee stability. 

A further major problem is that of the identification of instability, where the state of a net- 
work may remain over extended periods of time in each of a number of "quasi-equilibrium" 
distributions, some of which may correspond to highly inefficient performance. Instability 
is further closely linked to problems of phase transition in the probabilistic models of sta- 
tistical physics, and to the study of how phenomena such as congestion propagate through 
a network. At present results only exist for some very regular network topologies (see, for 
instance, Ramanan et al., 2002 and Luen et al., 2006). 

Questions related to those above concern the identification of fiuid limits, and in particular 
the problem of the uniqueness of their trajectories given initial conditions. It is notable 
that the uniqueness question has not yet been resolved even in the case of a general 
uncontrolled loss network, although it is known that here all trajectories do tend to the 
same fixed point, thus guaranteeing network stability. Further, while fixed points of fiuid 
limits identify quasi-equilibrium states of a network, detailed behaviour within such states, 
and the estimation of the time taken to pass between them, requires a more delicate 
analysis based on the study of diffusion limits. Here relatively little work has been done 
(see Fricker et al, 2003). 

Finally we mention that loss networks may be seen as a subclass of a more general class 
of stochastic models, with state space Z;^ for some R and fairly regular transition rates 
between neighbouring states. Notably their analysis has much in common with that of 



19 



processor-sharing networks, in which calls again have a simultaneous resource requirement. 
A unified treatment is still awaited. 
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